Coherent Destruction of Photon Emission from a Single Molecule Source: A 

Renormalization Group Approach 
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The photon emission from a single molecule driven simultaneously by a laser and a slow electric 
radio frequency (rf) field is studied. We use a non-Hermitian Hamiltonian approach which accounts 
for the radiative decay of a two level system modeling a single molecule source. We apply the renor- 
malization group method for differential equations to obtain long time solution of the corresponding 
Schrodinger equation, which allows us to calculate the average waiting time for the first photon 
emission. Then, we analyze the conditions for suppression and enhancement of photon emission in 
this dissipative two-level system. In particular we derive a transcendental equation, which yields the 
non-trivial rf field control parameters, for which enhancement and suppression of photon emission 
occurs. For finite values of radiative decay rate an abrupt transition from the molecule's localization 
in its ground state to derealization is found for certain critical values of the rf field parameters. 
Our results are shown to be in agreement with the available experiments [Ch. Brunei et al, Phys. 
Rev. Lett. 81, 2679 (1998)]. 



I. INTRODUCTION 

Pump-probe experimental setup, used and studied in 
optical spectroscopy for years (see, for example, [j) plays 
a significant role in fluorescence spectroscopy of single 
molecules. When a monochromatic off-resonant electro- 
magnetic field is added to already present near-resonant 
one, a single molecule exhibits a great variety of quan- 
tum optical effects. Just a few to be mentioned in this 
connection, are the light-shift of the zero-phonon line, 
Autler-Townes splitting [3] , multiphoton resonances @] 
and rf Rabi resonances jj. The current interest in this 
field, within the context of single molecule spectroscopy, 
stems from the ability to use the slow rf field to manip- 
ulate the single molecule, in such a way that a single 
photon or a few photons are emitted per rf period. Such 
single photon spectroscopy, is likely to be the most del- 
icate measurement one can perform on a molecule using 
optics. Furthermore, such an approach might be useful 
for quantum cryptology, and quantum computing. 

Possible frameworks for the analysis of existing low 
temperature optical experiments on fluorescence spectra 
of a single molecule 5] are provided by a well known two- 
level atom picture. Among these frameworks is a descrip- 
tion of photon emission by a single molecule via optical 
Bloch equations. This approach has a long history in 
quantum optics, and in studies of nuclear magnetic reso- 
nance and electron-spin resonance [g|. Quite frequently, 
however, one has to settle for the numerical solution to 
compare Bloch equations description of a single molecule 
with laboratory measurements. 

Another common analytical tool in this context is the 
quantum jump approach [7]. It treats the two- level sys- 
tem (TLS) emitting photons as an open or dissipative 
system where the time evolution of the system is gov- 
erned by the Schrodinger equation with an effective non- 
Hermitian Hamiltonian operator. The system evolves 
from its ground state \g) to a superposition of excited 
state |e) and state \g), and jumps back into state \g) 



emitting a photon. Then the process starts over again. 
It is possible to characterize such jumps statistically, and 
therefore to get the statistics of photon emission by single 
molecule source. 

The objective of this paper is to study the average 
waiting time of the emission event by a SM interacting 
with laser and off-resonant rf field. Using the renormal- 
ization group (RG) method for differential equations jj| 
we carry out the perturbative solution of time depen- 
dent Schrodinger equation with effective quantum jump 
Hamiltonian. Our goal is to obtain the conditions for 
enhancement and suppression of photon emission, as the 
rf field control parameters are changed. As discussed in 
Ref. 0|, this problem is related to the problem of co- 
herent destruction of tunneling (CDT) for a particle in a 
double well potential [lfl |. 

From CDT studies it is known that specific rf field pa- 
rameters may cause particle's localization in one of the 
wells |l(l |. when the transition probability between the 
two states of the system is made vanishingly small. Such 
a suppression or coherent destruction of tunneling |l7j 
corresponds to a suppression of emission in a SM 
However, the SM system is fundamentally different then 
an electron in a double well potential, since the SM emits 
photons. Therefore, we are left with the task of estab- 
lishing the influence of radiative life time of the molecule 
on the existing CDT picture. 

This is done in Section II, where we present the RG 
analysis for the equations governing the evolution of the 
dissipative TLS. In Section III we derive the criteria for 
localization of the system in its ground state for the cases 
of small and finite radiative decay rate and discuss the 
transition from localization to derealization as well as 
critical behavior of the spectrum. The conclusions are 
given in Section IV. A brief summary of our results was 
given in (Tll |. 

Finally, in connection to the fluorescence peak suppres- 
sion condition it is worth mentioning the well known re- 
sult from the theory of dressed states in quantum optics. 
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The intensity of the resonances of the atom dressed by the 
rf photons is given by Ji (0 (/ = 0, 1,2 . . . ), where the 
modulation index £ is equal to the inner product of sys- 
tems dipole moment and external field in units of driving 
frequency (see details below). Therefore, the weight 
of the resonance can be made infinitesimal, by tuning the 
modulation index to one of the zeroes of the correspond- 
ing Bessel function. For the connection between the the- 
ory of dressed states and perturbation theory used in the 
present paper see Ref. [12l |. Note that the I = case is 
precisely the CDT condition ^j. In this paper we find a 
new condition for the peak suppression which takes into 
consideration the radiative part of the Hamiltonian and 
modifies "zeros of Bessel function" condition, well known 
in quantum tunneling and quantum optics. 

II. RG ANALYSIS OF THE SCHRODINGER 
EQUATION FOR A SINGLE MOLECULE 

For a molecule with a ground electronic state \g) and 
an excited electronic state |e), in the presence of a laser 
and a slower rf fields, the Schrodinger equation has the 
form 

4 1* (*)> = #(*)!*(«)>, 



H (t) =V g cosu; rf t \g) (g\ + - (\g) (e\ + |e) (g\) 
+ ^V e cosidrft — + S^j |e) (e| . 

Here u> r f is rf field frequency, 5 is the laser detuning, 
V e .g = fi e g -E r f, Q = A*eg-E (the Rabi frequency), where 
/x e g are permanent dipole moments of the molecule in 
states |e) and |<?}, /x eg is the transition dipole, and E 
and E r f stand for the amplitudes of the laser and the 
rf fields respectively. The emission event is taken into 
account by adding a non-Hermitian part — {iT/2) |e) (e 
to the TLS Hamiltonian ;(j | . 

i^|*(*)) = (ff(*)-(»T/2)|e>(e|)|*(t)>, (1) 

where T is the rate of radiative decay, also referred to as 
damping coefficient. We are interested in the evolution 
of |* (*)) = * s (t) \g) + * e (t) \e), or, to be precise in the 
evolution of the survival probability 

p (*) = (*(*) !*(*)), 

i.e. the probability of no emission event to occur in the 
time interval between and t. Note, that in Eq. Q 
the TLS Hamiltonian H (t) is taken in the rotating wave 
approximation for the fast on-resonant laser field, but not 
for the slow off-resonant rf field. 

From the survival probability the average waiting time 
(t) for emission event is readily obtained, according 



to 

f°° f°° / \ 

(r)=J P Q (t)dt = J (\* g (t)\ 2 + \* e (t)\ 2 )dt. (2) 

This quantity yields the average waiting time for the first 
emission when (0)} = \g), and the phase of the rf field 
(Eq. QJ) is initially equal to zero. It can serve as an es- 
timate for the mean time between successive emissions, 
provided we can neglect the effect of the phase of the rf 
field. When the Rabi frequency is small enough, such an 
effect can be minimal, because the excited state popula- 
tion is proportional to fl 2 , when only laser field is present. 

In what follows, we discuss two cases: T/uj r f <C 1, 
n/uj r f ^ 1 and fl/uj r f <^ 1, T/w r f ~ 1. The first one, 
allows the comparison with fluorescence measurements 
carried out in Ref. and illustrates the localization 
condition of a single molecule in its ground state, induced 
by the rf field. In the second case the localization gets 
destroyed with the increase in dimensionless damping co- 
efficient T/uj r f, and we are able to follow the transition 
from localization to derealization. 

We should also note, that Eq. Q has time-dependent 
coefficients, and no closed form solution is known for a 
general set of laser and rf fields parameters. The simpli- 
fication introduced by choosing a small Rabi frequency 
in Eq. (fl/to r f <C 1) leads to the possibility of per- 
turbative calculation of (t)}. 

A. RG analysis in the limit of small radiative 
decay rate 

First of all we select the time scale of our problem and 
nondimensionalize all its parameters. This can be done 
as follows: t i— > oj r ft, V e>g i— > V e ,g/w r f, 6 5/u r f, T i— > 
T/u> r f, Q i ► Q,/bj r f. By referring to the Rabi frequency 
as being small we mean, that it is small compared to the 
frequency of the rf field. We can further simplify our 
analysis of Eq. by assuming the system dissipation 
to be small as well. More explicitly, T i— > eT, fi i— » efi 
(e«l)- 

Next step is the change of variables accord- 
ing to: f g (i) = Cg(t)exp{-iV g smt}, <J e (t) = 
c e (t) exp {—W e sin t — iSi}, to prepare for the i nteg ration 
in Eq. Then, we make use of the identity 

oo 

exp{izsint}= V" e ikt J k {z), (3) 

k— — oo 

and the original system (Eq. QJ) reduces to the following 
set of coupled ordinary differential equations: 

|c s = -z £Ce ^ £ e-"j fc (£), 

k— — OG 

±C= -iec g ^ £ e ikt+lSt J k ({J - 6 l Ce (t) , (4) 

A;— — og 
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where we introduced the modulation index £, — V e — V g . 

We now turn to the perturbation analysis, with the 
objective to get an asymptotic solution of Eqs. Q) which 
would be valid globally in time. To do that we apply 
the RG method introduced in B. Unlike conventional 
singular perturbation methods [1J|, the RG procedure 
starts with the construction of the naive expansion. It 
takes the form of a power series in the small parameter 
e: 

C (t) = C (°) (t) + ec« (t) + e 2 c^ (t) + O (e 3 ) . (5) 

The naive solution is then obtained by substituting C (t) 
into Eq. J3J and solving for the gauge functions (t) 
at each order e n . It is assumed, that the gauge functions 
are 0(1) for all t. However, this is almost always not 
the case. Usually, for large enough times one or several 
c'") (t) become greater than their predecessors v-" 1 ' (t) 
(m < n). Then, the naive expansion becomes nonuni- 
form and cannot be used as an approximate solution. 
The terms causing nonuniformity are called secular terms 
[l^j . The RG procedure regularizes the naive expansion 
by identifying the secular terms and eliminating them. 



Although, there exist a bulk of literature on RG method 
and its applications 0], we outline its key steps and 
details in Appendix A. 

In the naive expansion (Eq. (2J), which we now sub- 
stitute into Eq. J3J, all the gauge functions are two- 
component vectors c 1 -™' = fcg"\ci"^ . At O(e ), the 

solution of Eq. Q is a constant vector c' ' = (A,B) T . 
One can see, that at O (e) this leads to the linearly grow- 
ing terms (secular terms) at integer values of detuning: 
5 = 1. Thus, to be able to treat other values of detuning 
simultaneously with the integer ones, we set 8 = — I + ea, 
where a — O (1) is the "detuning of the detuning" and 
proceed. The solution we are about to obtain describes 
the vicinity of fluorescence peak (cf. Ref. Q). The entire 
spectrum can be approximated as a sum of such peaks, 
located at different integer values of detuning, provided 
that the peaks are not too wide. Since we consider the 
small dissipation and small Rabi frequency case, this as- 
sumption is secure. 

We introduce a parameter a into Eq. <QJ to get 



dt 



£1 s.itecr 



k^l 







(IV 1 ikt+iU 



■MO 



-ikt- 



2 / , 

k^l 



a 



(6) 



with C = (c g ,c e ) T , and C (0) = (1,0) T . From Eq. © 
the result to O (e) reads 

A-^J-t (0 (t-t ) 



C = 



B 



J_ ; (0 (i - to) - (e£ + io) B(t- t Q ) 



where to is an arbitrary time scale used to regularize the 
naive expansion (See Appendix A). In the O (e) solu- 
tion we dropped all the terms other than secular, due 
to the irrelevance of such terms for the next steps of the 
procedure. The RG equations are obtained by imposing 
dC/dto = at to = t- We have: 



d_ 

dt 



Br 







%J-i (0 



.4 



where subscript U R" indicates, that we are solving for 
the renormalized values of A and B. Finally, solving Eq. 
with the initial conditions C (0) = (1,0) T we obtain 
the following renormalized zero order solution C (t) = 

(A R (t) , Br (t)f + 0(e): 

Ai exp {e\2t} — X'2 exp {eAit} 



Ar (t) 



Ai — A2 



(8) 



Br (t) 



Ai-A- 



(exp{eA 2 t} - exp {eAii}) , (9) 



Ar, 



- - — ± -Jv 2 - 16A2 - 4a 2 - AiaT. 
4 2 2 V 1 



Here, by "±", we specified two different square roots, and 



We can now set e = 1, as we do not need it anymore. 

After the integration in Eq. (Jgjl, it can be shown, that 
(t} ; , the average waiting time, when the detuning is close 
to the integer value I, is given by 



- 8A 2 + 4ct 2 
IfA 2 ' 



Then, for arbitrary laser detuning S, the average waiting 
time, can be approximated as 



v-^r 2 + 8A 2 + 4a 2 



1 



4rA? 



(10) 



The equality sign in Eq. IllOfl indicates the approximate 
nature of this formula, which we mentioned above. The 
first thing that catches attention is the presence of Bessel 
functions in the denominators of each term in a sum. 
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This opens the possibility to suppress the fluorescence 
peak (to infinitely extend waiting time), and to localize 
the single molecule in its ground state. For integer values 
of the detuning S — I (in units of w r /) we can always 
select parameters of the rf field, i.e. modulation index £, 
so that J s (0 = 0. 

In Fig. Qk, we compare the RG result (Eq. ifTfljl ) 
with the numerical solution of Bloch equations, previ- 
ously done in Ref. Q in order to match it with their ex- 
perimental results. We repeat numerical simulations of 
Ref. i for the same set of system parameters and calcu- 
late the fluorescence spectra. The prediction of Eq. (1101) 
is in a good agreement with the numerical solution for 
the three smallest values of Rabi frequency (Fig. Q . The 
slight discrepancy observed for the highest value of Rabi 
frequency (top graph in Fig. QJi) is due to the fact that 
perturbation parameter e « 0.46 in this case is rather 
close to 1. In particular, the RG solution (to this order) 
does not capture the shifts in the fluorescence peaks po- 
sition, found in experiments for CI > 3.2r U In Fig. 
Q]d we demonstrate the outcome of changing the value of 
modulation index £ from 1.14 of Ref. [5] to the first root 
of zeroth order Bessel function (£ ~ 2.4048). We can see 
that the criterion for destruction of emission (Js (£) = 0) 
is supported by numerical solution of Bloch equations, 
as the central peak (5 = 0) gets suppressed for the latter 
value of £. 

Now we further explore the emission suppression case. 
Keeping the same parameters as in Ref. Jj| we change 
value of £ to make the side peaks disappear. The results 
are presented in Fig. We also observe, that for suffi- 



Proceeding with the naive solution, we observe that the 
secular terms do not appear in this expansion at 0(e), 
as they did in the case of weak dissipation. 

At 0(1), the solution is a constant vector cS°^ = 
(A,B) T . Substitution of this result into the O (e) equa- 
tions yields: 

k— — oo 



c« = iCIA ± J k (0 (e»/*H* _ i) 

k— — oo 



ciently large £, i.e. large number of the Bessel function 
root, we can suppress more than two side peaks, cor- 
responding to positive and negative integer detunings. 
Taking £ to be the 8th root of J\ (£) we eliminate the 
third peaks, while taking it to be the 8th root of J 2 (0 
we eliminate the central peak (S = 0) as well. This can be 
explained by the behavior of zeroes of Bessel functions, 
which fall densely on the real axis, as the argument in- 
creases. For the large modulation index one expects the 
entire spectrum to be suppressed (l^ . 



B. Average waiting time for the case of finite decay 

rate 

The evolution of survival probability alters signifi- 
cantly when the dissipation coefficient T is larger, or of 
the same order as the driving frequency u) r f . This be- 
comes clear from the structure of the secular terms as we 
take on the long time asymptotic solution for this case. 
We still consider Rabi frequency to be small (fl *—> eO, 
e<l), but keep T = O (1). First of all, we focus on the 
case of zero detuning (S = 0) . The Schrodinger equation 
then reads 

il#=( v * co , st en / 2 , U (id 

dt \ V e cost-iT/2) ** {LL) 

The transformation to "amplitude" variables is almost 
the same as before: Vg(t) = c g (t) e~ iV ° sint , * e (t) = 
c e (t) e~ lV " smt - rt / 2 . Once again, we face a set of coupled 
linear ODEs: 



(12) 

I 

which, in turn, leads to the following O (e 2 ) equations: 

k,l— — 00 

( vt/2+at _ 1 A r - 2d 

k,l— — 00 

x d e -rt/2-ut\ r-2il 
Further integration produces secular terms (for k = I) 

9 2 ^ T2+4F + ' 

A.: — — OO 



dt 



n 

-i — e 
2 





Jkt+1 



Eoo 
fe=- 



Jk(0 



c. 



5 



(2) «7_T J fc (Q 



k=—oo 



where "NST" stands for nonsecular contribution to O (e 2 ) 

solutions. The RG equations (to order 0(e 3 )) for the f oo 1 

regularized "constants" A and B are: ^ (t) — B (0) ex < 6 — ^ ^ fc (0 I /.« 



= ~ A n- r - E 



dt " 2 ^ r 2 + 4fc 2 

fc— — oo 



r 2 + 4fc 2 

k=—oo 



d_Bn = ~ J 2 (Q 

dt R 2 ^ r 2 +4/c 2 ' 

&J — — oo 

Hence, to O (e 3 ) A# and satisfy 



At time ( = we assume that the molecule is in the 
ground state: c e (0) = 0, c g (0) = 1. Applying these ini- 
tial conditions and replacing A and B with their renor- 
malized values (Eqs. \1?>1 ) in we obtain O (e 2 ) per- 
turbation results for amplitudes c e and c g . Then, we 
switch back to the original variables ^> e , g (t) and arrive 
at 



^ _ e -iV g sin t 



exp ■ 



e 2 ^ 2 r ~ 4(0 



2 ^ r 2 + 4fc 2 

k— — oo 



l-e 2 Sl 2 £ 



■MO 4(20 fi(l-e^- fe )*) 2(1 



k^l — — oo 



T + 2U 



I - k 



T + 2ik 



(14) 



* P = -ieSle 



—iV e sin £— 4-i 



exp 



E 



oo 

E 



Tt/2+ikt 



Jt (0 

2 ^ r 2 + 4fc 2 

k— — oo 

r - 2ifc 
r 2 + 4/c 2 ' 



(15) 



Then, using Eq. (2J we calculate the mean waiting time 
(r); the O (e 2 ) result reads 



(t)- 1 = rn 2 



E 

/c=— oc 



j2 (0 

2 +4fe 2 



Q2\ 
fa J 



(16) 



The corresponding result for the case of nonzero detun- 
ing can be calculated in a similar fashion. Formally, this 
case is different from the one we just considered by ad- 
dition of imaginary part to the decay rate: r i— > T + 2iS. 
Therefore, the general formula for the case of arbitrary 
laser detuning can be written as 



(rr x = ro 2 £ 



4 (2?) 



The possibility of locking the system in state \g) for a 
long time, is still present. As one can see, the presence of 



finite damping has only modified the corresponding con- 
dition, and made (r) finite at all values of £ in contrast 
to the results of Eq. IIOII . In other words, the average 
waiting time can be maximized for certain values of mod- 
ulation index, but the height of the peaks at the maxima 
positions is now finite even at the leading order in small 
Rabi frequency. In Fig. El present the comparison of this 
prediction with the results of numerical solution of Eq. 
(QJ. Predictions for all three values of T are adequate. 

From Fig. El we observe, that when the dimensionless 
r is not too large, the maxima in the average waiting 
time occur at zeroes of </<5(£), m accord with the re- 
sults obtained in small dissipation limit. If the decay 
rate is increased, one expects the peaks to get broaden 
and shifted. This result does follow from the comparison 
of T = 0.5 and T = 1.0 cases (see Fig. EJ). Meanwhile, 
the third curve in Fig. corresponding to the case of 
r = 3.0, points to a less trivial feature in the peaks be- 
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havior. Considering the presence of well defined maxima 
in (t) as a sign of localization (in ground state) effect, 
we can refer to their disappearance as a derealization 
effect. 

The detailed study of such a transition is the purpose 
of the next Section. 



tion of corresponding condition in absence of dissipation). 
When dissipation becomes noticeable, the peaks of (r) 
shift with respect to the £„. Assuming that these shifts 
are still small, one can estimate them perturbatively from 
Eq. iJTEJ. 



III. TRANSITION FROM LOCALIZED TO 
DELOCALIZED REGIME 

As long as the decay rate does not exceed certain 
threshold value, the average waiting time (r), plotted 
as a function of £, appears as a sequence of alternating 
minima and maxima (see Fig. EJ. Those extremums 
correspond to the enhancement and reduction in photon 
emission rate, respectively. Setting d (r) /c?£ = and dif- 
ferentiating Eq. itlfi)) we obtain the following conditions 
for the extremums: 



■MO 



Ji(0 



E 



J k (0 ( Jk-i (0 - Jk+i (0) 
r 2 + 4/c 2 



(18) 



G(r,o 



Condition for the emission suppression in the 
small shifts limit 



At the position of a maximum, we have d (r) /<i£ = 0. 
Upon differentiation in Eq. itlfi)) , we use the sum rule 0] 

Er = -oo 4 (0 = ^EZ-ooJi (0 = i- J 2 o (0, which im- 

plies 



2][>(0 



k=i 

and find, that 



dJk (Q 



-Jo (0 



dJp (Q 

da • 



d(r) 
d£ 



E 

k=i 



■4(0 



r 2 + 4fc 2 



Thus, for zero laser detuning, if T is infinitesimal, the 
nth maxima would be located at £„, the nth zero of the 
Bessel function J (£) (see also Appendix A for the deriva- 



Then we use an 

(Jfe-i (0 - Jk+i (0) /2 



identity dJk (0 /d£, 
and addition theorem [H 



J k (U±V) = J2hL-oc J kTl ( u ) J i ( v ) ( M < M) t0 S et 



J 



d£ ^ 



Eoo 
l=- 



Jk-l (£n) Jl (v u )Ya= 



X 

/ — — oo 



{Jk-l-l (Cn) - Jfc-i+1 (£«)} Jl {v n ) 2 



r 2 + 4fc 2 



where v n is the shift of the nth peak. Next we make use 
of the small v n asymptotics of Bessel functions: Jo { v n) ~ 
1 - ^ 2 /4, Ji (v n ) ~ {v n /2) 1 (I > 1). Keeping the zeroth 
and the first order terms, and solving for v n we arrive at 

_ Dfcll r 2 +4k 2 J k - Jk+l (£n)} 



with 



(19) 



p k = {Jk-i - Jfe+i (e,)} 2 , 

Qfc = Jfc (£„) {4-2 (Cn) - 2J fe (£„) + J fc+ 2 (£„)} 

For n = 1 through 4, we plot the relative shifts v n /£n> 
given by Eq. iflfll) , as a function of damping coefficient T 
in Fig. 31 and compare them with numerical solution of 
the transcendental Eq. 1181 . 

It is clear that estimates given by Eq. Q19JI predict the 
shifts, which are less than 10%, but fail to reproduce the 
larger shifts. Further analysis of extremums behavior has 
to be done by solving Eq. II18|I numerically. 



Condition for the emission suppression for 
arbitrary shifts 



In Fig. 0we plot G(T,£) and Jo (0 versus £ to illus- 
trate the solution of the transcendental Eq. itTsl) for two 
values of r. When T = 1.0 the corresponding curves cross 
close to £ n (Fig. However, as T is increased to 2.5, 
Eq. l|T8jl does not have any solutions in the vicinity of £i 
and £2 and the first two maxima of (r) disappear (Fig. 
0a). We also notice, that for this value of the decay rate, 
the first two minima of (r) , located next to the zeroes of 
Ji (£) (Eq. d} ) f° r smaller T, disappear as well. 

Apparently, maxima and minima of the waiting time 
(r) do not vanish gradually due to the peaks broadening. 
Instead, there is a critical value of T (T/u) r f in original 
units) for each extremum. Therefore, with an increase in 
dimensionless decay rate T/uj r f a single molecule under- 
goes the transition from localized-like behavior to delo- 
calized behavior, for which we can no longer speak about 
maximizing average waiting time, or about the localiza- 
tion of a molecule in its ground state. The transition 
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manifests itself as a motion of the neighboring pairs of 
maxima and minima towards each other with subsequent 
annihilation at the same values of critical damping T cr 
(see Fig. EJ. The shifts, given by numerical solution of 
Eq. I|18J1 are plotted versus radiative decay rate T in Fig. 
H3 They are in accord with the shifts obtained from the 
numerical solution of the Schrodinger equation for the 
first two peaks of average waiting time. Slopes of both 
curves in Fig. Ujjstart to diverge, as the dissipation grows, 
in anticipation of critical points. For example, r cr ~ 1.7 
for the first peak. 

At a critical point, the slopes of Jo (£) and G(r,£) in 
Fig. coincide. For the dissipation approaching its crit- 
ical value we can use this fact and determine the asymp- 
totic behavior of extremums. Expansion in Taylor series 
of the Eq. I|18H in the vicinity of critical point (£ cr , r cr ), 
produces 

d 2 Jo fcr) (t _t \2 dG (r cr , £ er ) 

de 14 Wj ~ dv [ cr) - 

Here we retained only the leading terms in both sides 
of the equation. Thus, the position of extremum scales 
with the decay rate as £ cr — £ oc (r cr — T) 13 , with f3 = 1/2. 
Further illustration to the transition from localization to 
derealization is provided in Fig. Q where we display the 
graphical solution of Eq. l|T8jl in (£, F)-plane. 

IV. DISCUSSION AND CONCLUSIONS 

We have shown, that a dissipative two-level system, in- 
teracting with continuous laser field can be localized in its 
ground state by appropriate choice of parameters of ad- 
ditional driving field. To demonstrate this we considered 
a single molecule governed by a non-Hermitian Hamilto- 
nian in rotating wave approximation for the laser field. 
The added slow rf field interacting with a molecule via 
dipole moments of its ground and excited states makes 
the TLS Hamiltonian time dependent. Our main quan- 
tity of interest, the average waiting time for the first emis- 
sion requires knowledge of the survival probability, which, 
in turn, demands the solution of Schrodinger equation. 

Generally, this goal can be achieved only by pertur- 
bation analysis in the case when some of the governing 
parameters are small. For this reason, throughout the 



paper we assumed that the Rabi frequency Q of a sin- 
gle molecule is small compared to the frequency of the 
driving rf field w r f- On the other hand, to provide a 
correct estimate for the average waiting time, the result- 
ing survival probability has to be given by a uniformly 
convergent series. In this paper we overcame this diffi- 
culty by applying the RG procedure formulated in Ref. 
0. For the radiative decay rate of the excited state F, 
we considered two limiting situations: T/uj r f <c 1 and 
T/uj r f ~ 1. The RG calculation for the first case gener- 
ates two important results. 

(i) The fluorescence spectra of single molecule, calcu- 
lated as a reciprocal average waiting time reproduces ex- 
perimental results of Ref. jj| . 

(ii) There is a condition for suppression of fluorescence 
peaks, which resembles the criteria of coherent destruc- 
tion of tunneling ^3] of electron in a double well. 

Similar procedure in the second case also produces the 
correct estimates for the average waiting time (r) , which 
was verified by numerical solution of Schrodinger equa- 
tion for a single molecule. However, the localization of 
a molecule in its ground state is achieved for the values 
of modulation index different from those in the small T 
limit (the well known condition Jo (£) = does not ap- 
ply). Moreover, at a certain critical value of dissipation 
coefficient, the neighboring minima and maxima of (t) 
annihilate, and the localization no longer persists. 

In connection with the transition from the localized 
behavior to the delocalized one, we were able to estimate 
the shifts in the positions of extremums of average wait- 
ing time for a given dissipation, and to determine the 
critical values of T/uj r f for a given extremum. 

Finally, we remark, that similar analysis can be carried 
out for the system of optical Bloch equations, which rep- 
resents a standard tool in single molecule spectroscopy, 
nuclear magnetic resonance and electronic spin reso- 
nance. 
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Appendix A: RENORMALIZATION GROUP 
METHOD FOR DIFFERENTIAL EQUATIONS 

Here we illustrate basic steps of the RG method Q on 
two examples. The first one is a trivial linear differential 
equation 

dtip + iuJTp = eif), 

with e <C 1. Although the exact solution is known, we 
wish to recover it perturbatively. Below, and everywhere 
in the text, the regularization of the naive solution is 
achieved by analyzing the integration constants enter- 
ing the solution. These "constants" play the role of bare 
quantities of quantum field theory and statistical me- 
chanics. They are treated as functions of time and are 
modified (renormalized) in order to remove secular terms. 



We start by looking for a solution given as a power 
series in small parameter e: 

i/> = + #1 + e 2 ^2 + O (e 3 ) , (Al) 

where i/'fe are O (1) quantities. At lowest order, we have 

dt^o + = =h ipo = A exp {— iuit} . 

We substitute this result into the equation for the first 
correction to get 

dttpi + iojipi = A exp {— itot} , 

which leads to 

tpi = At exp {— icut} , 

as we only need a special solution of the inhomogeneous 
equation (the solution to the homogeneous one is ab- 
sorbed into the zeroth order perturbation result i/'o)- The 
first correction grows linearly, and after 1/e amount of 
time is elapsed Vi becomes comparable to supposedly 
dominating term -00- Thus, 

tfj = A exp {-iut} + etA exp {-iujt} + O (e 2 ) . (A2) 

does not represent uniformly valid result. Nevertheless, 
according to the RG formalism, this expansion gives us all 
we need for the construction of O (e 2 ) long time asymp- 
totic expansion. 

We now assume that A is not a constant, but a function 
of some arbitrary time scale to- The renormalization con- 
sists of two stages: first, we remove divergency to O (e 2 ) 
by the following manipulation: 

-0 = A (1 + eto) exp {— iuit} + e (t — to) A exp {— icot} 
= A R (t ) exp {-iuit} + e (t - to) A R (to) exp {-icut} , 

(A3) 

literally "renormalizing" A, and introducing new variable 
A R instead. The latter expansion, is uniform in time, as 
r can be chosen arbitrarily large. However, -0 should not 
be a function of some artificial time scale. Therefore, at 
the second stage, we require 

J-V (to =t)=0 (A4) 

Oto 

Next, we eliminate to by setting it equal to t (since to 
is arbitrary, it is at our disposal). As a result of Eq. 
IIA4j) . we obtained differential equation (called the RG 
equation) governing our former constant A. It reads: 

±A R (t) = eA R (t) + 0(e 2 ). 

Finally, solving the RG equation, substituting the result 
into the expansion <|A3|) . and setting r = t we arrive at 
the long time asymptotic solution uniformly valid to the 
next order in e: 

tp = A (0) exp {et -iut + O (e 2 t) } + O (e 2 ) . 
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For this simple problem it can be further shown, by 
repetition of the above procedure at all orders, that 
if) = A (0) exp {et — iujt] is actually a global solution. 

A more complicated example of Schrodinger equation 
for the electron in a double well potential ^3] can be 
also treated using RG method. Here, our presentation is 
slightly different from solution of the this problem given 
in Ref. ^3], but it leads to the same results and illus- 
trates all the main steps of the derivations in Section II. 

The TLS Hamiltonian in this case, is: 

H = V g cosut \g) (g\+V e cosut \e) (e|+- (\g) (e| + |e) (g\) . 



c e = B R - ieA R ^ (t - r) J (0 + NST + O (e 2 ) , 

where NST stands for the nonsecular terms. As the next 
step, we obtain following RG equations: 

^ = -^W0+O(e*), 
= --A R J (0 + O(e 2 ). 

The latter system is readily solved, and with the initial 
conditions c g (0) = 1, c e (0) = produces: 



With this Hamiltonian, the Schrodinger equation for the 
TLS, which is initially in its \g) state results in the fol- 
lowing system of coupled ordinary differential equations 
for amplitudes c e and c g 



( 



dC_ 
dt 



iQe 
~2~ 



E 



-ikujt 



\ 



\ k— — oo 



■MO 



c 



(A5) 

where C — (c g ,c e ) , and C (0) = (1,0) and we used 
Eq. Q. 

At O (e°), the solution is simply C (0) = (A, B) T , with 
both A and B being constants. Naively, we would apply 
the initial conditions at this stage. Instead, we substitute 
this result into the equations for next order correction, 
and obtain 



.(i) 



-i^B{J (0t 



E 

k——-oo,k^0 



Jk (0 (e- 



ih-l 



+ O (e 2 ) , (A6) 



c e = sm <j —J (0 n + e- cos i — J (£) t 



x E 

k— — oo,k^0 



1-e 



OU 2 



(A7) 



Note, that coherent destruction of tunneling condition 
Jo (2£) = [10| follows immediately from these results 
(see also Ref. 12]). Continuing the procedure to order 
O (e 2 ), we obtain 

dA R _ 2 n 2 ietl 



dB R ieVt . , . 9 fl 2 . . . 

= ~— ^.Jo (0 - e 2 T FB fl , (A8) 



c p=-i^A{J (0t 



E 

A;— — oo, k^0 



Jk (0 (i - e 



fcn 



The renormalized asymptotic expansions are 

c g = A R - ieB R ^ (t - r) J (0 + AST + O (e 2 



with 



* = -*■*>«) E *f-< E -' ; 

fe=— oo,fe^Q k— — oo,k^0 



4 (0 



Even without solving Eq. l!A8|l . we observe, that lo- 
calization condition holds to order e 3 . For the recently 
found third order corrections see Ref. 0] and references 
therein. 
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Figure 1: Comparison of the predictions for {t} -1 (5) given 
by Eq. 11011 (solid lines) with the excited state population 
calculated from Bloch equations (empty triangles) (Ref. 0|) 
for various values of Rabi frequencies (from bottom to the 
top): fi = 0.29; 0.9; 3.2 in units of Y. Parameters are taken 
from experiment of Ref. |J]: Lu rf /2-7v = 140 MHz, (a) £ = 
1.14; (b) £ is the first root of Jo (£) = 0, F/2tv = 20 MHz. 
The results are arbitrarily shifted vertically for transparency, 
therefore, units on vertical axis are arbitrary. Our predictions 
in (a) are in accord with experiments p|. In (b) we observe 
the suppression of the fluorescence for 5 = 0. 
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Figure 2: Suppression of the fluorescence peaks for different 
values of modulation index. Here = 0.9F, r/27r = 20 MHz 
(parameters of Ref. Q). (a) £ is the first root of Ji (f) = 0; 
(b) £ is the first root of J2 (£) = 0; (c) £ is the eighth root of 
Ji (£) = 0; (d) £ is the eighth root of J 2 (£) = 0. 
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Figure 3: Scaled mean time between emission T (r) is plotted 
as a function of modulation index £ for the case of zero detun- 
ing. The solid line represents RG prediction of Eq. llfit . dots 
correspond to the numerical solution of Eq. Grid- lines 

indicate the positions of zeroes of Jo (£). The Rabi frequency 
is = 0.1. The radiative decay rates are F = 0.5; 1.5; 3.0 
(from top to bottom). We notice the shifts in maximum of 
(r) relative the zeroes of Jo (£). 





Figure 5: Illustration to the transcendental Eq. I|18jl . Bessel 
function Jo (solid line) plotted as a function of £ together 
with the right hand side of Eq. 11811 . G(T, £). Radiative 
decay rates are (a) T = 1.0, (b) T = 2.5. Grid-lines indicate 
zeroes of Jo (£). In localization case (a), the solutions of Eq. 
11811 are slightly shifted with respect to £ n . In derealization 
case (b), these solutions do not exist since the curves do not 
cross. 
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Figure 6: The shifts (£ - £„) /£„ found from Eq. (O (solid 
lines) are compared to those found from numerical solution of 
Eq. £Q| (boxes). They are plotted as a function of F (T/ui r f 
in original units) for the peaks close to the first (upper curve) 
and the second (lower curve) zeroes of Jo (the first two peaks 
in Fig. 01. Critical points are marked by stars. They indicate 
values of T /u> r f above which the transition to derealization 
occurs. 
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Figure 7: Graphical solution to Eq. 1181 represented as a 
solid curve in (£, F)-plane. The solution represents the value 
of control parameters which maximize and minimize (t) . The 
left edge of each "tooth" corresponds to the maxima (inter- 
secting the £-axis at the zeroes of Jq (£): 2.40483, 5.52008, 
8.65373, etc.); the right edge corresponds to the minima. The 
critical points (we can see three of them) are identified as 
"teeth" tips (cf. Fig. 0. 



